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We investigate membrane-mediated interactions between transmembrane proteins using coarse- 
grained models. We compare the effective potential of mean force (PMF) between two proteins, 
which are always aligned parallel to the z-axis of the simulation box, with those PMFs obtained for 
proteins with fluctuating orientations. The PMFs are dominated by an oscillatory packing-driven 
contribution and a smooth attractive hydrophobic mismatch contribution, which vanishes if the 
hydrophobic length of the protein matches the thickness of the membrane. If protein orientations are 
allowed to fluctuate, the oscillations are greatly reduced compared to proteins with fixed orientation. 
Furthermore, the hydrophobic mismatch interaction has a smaller range. Finally, we compare the 
two-dimensional thickness profiles around two proteins with the predictions from the elastic theory 
of two coupled monolayers, and find them to be in very good agreement. 



I. INTRODUCTION 

Biological membranes are essential components of all 
living organisms. Even though their basic building block 
is a fluid lipid bilayer, their properties and functional- 
ity crucially depend on the membrane proteins. Pro- 
teins function as catalysts, transport and store other 
molecules, provide mechanical support and immune pro- 
tection, generate movement, transmit nerve impulses, 
and control growth and differentiation [T] . 

Unfortunately, the components of real biomembranes 
are too diverse and complex to obtain detailed and un- 
ambiguous information about interactions of transmem- 
brane segments with a lipid bilayer [2]. Moreover, struc- 
tural perturbations or transformations of the lipid bi- 
layer are among the most difficult processes to probe ex- 
perimentally -3i. Thus, complementary theoretical ap- 
proaches and computer simulations of membrane sys- 
tems of well-defined compositions are necessary to elu- 
cidate the role of the lipid bilayer in processes like pro- 
tein aggregation and function. In this context, the use of 
coarse-grained models has become popular [1HI|. Even 
with today's computing resources, atomistic modeling of 
multicomponent lipid bilayers on length scales of several 
nanometer still remains a challenge. Moreover, coarse- 
grained model simulations give insight into generic prop- 
erties and mechanisms in lipid membranes, which are dif- 
ficult to extract from fully atomistic simulations of spe- 
cific membranes. 

In this paper we focus on lipid-mediated interactions 
between proteins or, more generally, membrane inclu- 
sions. We present simulational and theoretical results on 
fluid membranes containing model transmembrane pro- 
teins, i.e. proteins that span through the membrane. 
Many theoretical models have been developed to get a 
better understanding of the membrane-mediated inter- 
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actions between such proteins [51-[T5]. In these models, 
the inclusions were modeled as straight cylinders, which 
are aligned with the bilayer normal. Only recently, com- 
putational studies with coarse-grained models have been 
performed on the same problem, where both proteins in 
single-component bilayers P^HTB] and the effect of choles- 
terol on protein-lipid and lipid-mediated protein-protein 
interactions have been studied ^9j. From these stud- 
ies, the following general picture has emerged: The lipid- 
mediated protein interactions can be divided in short- 
range and long-range contributions. The short-range 
contributions depend on the local structure of the lipid 
layer and reflect packing and layering effects. The long- 
range interactions result from the elastic distortion of the 
membrane due to the insertion of the proteins and can be 
tuned, by tuning the length of the hydrophobic section 
of the protein ("hydrophobic mismatch" interaction). If 
the hydrophobic length of the proteins matches the thick- 
ness of the membrane, they will vanish, otherwise they 
tend to be attractive. In the literature, a number of other 
mechanisms that would induce long-range lipid-mediated 
protein interactions have been discussed (20- 22] . but for 
straight cylindrical proteins, the hydrophobic mismatch 
interaction seems be dominant. 

The work presented here is based on a coarse-grained 
lipid model [23l [24] , which reproduces the main phases 
dSHlZj and elastic properties [H HT] of DPPC bilay- 
ers. In previous papers, this model was used to study 
the lipid-mediated interactions between infinitely strong 
straight cylinders with fixed orientation along the bilayer 
normal in stressfree (TS] and stressed membranes |27j . 
The rigid restriction was motivated by the huge amount 
of related theoretical work (see above). In reality, how- 
ever, proteins may tilt in the membrane, their orientation 
may fluctuate, and this affects the lipid-mediated protein 
interactions. 

The present paper focuses on this effect. We study 
the lipid-mediated protein interactions, using two vari- 
ants of a coarse-grained protein model. The flrst vari- 
ant is the infinitely long straight cylinder studied earlier 
PEl 127] . and the second variant consists of a cylinder 
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of finite length with full freedom to tilt. This allows to 
assess in detail the effect of orientation fluctuations on 
the protein-protein interactions. In addition, we charac- 
terize the profiles of membranes containing two proteins 
and compare them with the theoretical prediction of an 
elastic theory. 

Our paper is organized as follows: In the next sec- 
tion, we introduce the simulation models and method 
and briefly comment on the elastic theory with which the 
simulation data are compared. The simulation results are 
presented in Section III. We summarize and conclude in 
Section IV. 

II. MODELS AND METHODS 

In the following a brief summary of the lipid model 
and the protein models used in this work will be given. 
In our study, we vary both the hydrophobic mismatch 
of the proteins (i.e., the length of the hydrophobic sec- 
tion), and the hydrophobic strength. Since hydrophobic 
interaction does not arise from the binding of nonpolar 
molecules to each other but from preventing polar sol- 
vent molecules from achieving optimal hydrogen bind- 
ing, the strength of the hydrophobic interaction depends 
on the relative polarity of both the solute and the sol- 
vent 28 . Experimentally the hydrophobicity of proteins 
can therefore be tuned by changing amino acid residues 
with different hydrophobicity of the protein [5]. Alanine, 
e.g., is less hydrophobic than leucine [29j. Alternatively, 
changing the pH of the solvent and thereby making pri- 
orly uncharged side-chains charged, will also affect the 
hydrophobic interaction of the lipid bilayer and the pro- 
teins [50] . 

A. Lipid Model 

Lipid molecules are represented by chains consisting of 
one head bead and six tail beads, and there are additional 
solvent beads |24j. Within the lipid chain, adjacent beads 
at a distance r interact via a finite extensible non-linear 
elastic (FENE) potential 

^fene(7') = -;^eFENE(Armax)^ln 1 -r-^ ) , 

(1) 

with the spring constant efene, the equilibrium bond 
length To, and the cutoff Armax- The angles 9 between 
subsequent bonds in the lipid are subject to a stiffness 
potential 

1^3A(e) =eBA(l + cos(?), (2) 

with the stiffness parameter eba- Beads of type i and j 
which are not direct next neighbors in a chain interact 
via a truncated and shifted Lennard- Jones potential, 

Vu{rla,,) - i ^y-^- j " ^-^'J' ^ < ^-^^^ (3) 
I otherwise. 



The offset Vc^ij is chosen such that Vi^j{r/(Tij) is contin- 
uous at the cutoff Tc^ij. The parameter aij — {ai + aj)/2 
is the arithmetic mean of the diameters ai of the interac- 
tion partners, and rc,ij = 1 Uij for all partners {ij) except 
(tt) and (ss): r^tt = '^o'tt and rc^^s — 0. Hence tail beads 
attract one another, all other interactions are repulsive, 
and solvent beads do not interact at all with each other. 

We use the model parameters j24^ ah = 1-1 cr*, t'o = 
O.Tat, Armax = 0.2crt, epENE = WOe/a'^, and sba = 
4.7 e. At the pressure P = 2.0 e/a^, the model repro- 
duces the main phases of phospholipids, i.e., a high- 
temperature fluid La phase at temperature ksT > 
ksTm ~ 1-2 e and a low-temperature tilted gel {L^/) 
with an intermediate modulated ripple {Pp') phase |25j . 
The energy and length scales can be mapped to SI- 
units [18] by matching the bilayer thickness or, alter- 
natively, the area per lipid and the temperature of the 
main transition to those of DPPC, giving 1 dt 6 A and 
1 e '-^ 0.36 X 10~^" J. The elastic properties of the mem- 
branes in the fluid state were then also found to be com- 
parable to those of DPPC membranes )l8j . 




FIG. 1: Snapshots of the protein models considered in this 
work. Left: Infinite cylinder with fixed orientation. Right: 
Finite tiltable cylinder. L denotes the total length of the 
hydrophobic section. The dark shaded band in the middle 
indicates the region where the hydrophobic attractive inter- 
action with lipid tails reaches full strength, the two adjacend 
light shaded bands show the the region where it smoothly 
drops to zero. The finite cylinder (right) is effectively capped 
by repulsive hemispheres (not shown). 



B. Protein Models 

1. Cylinder of infinite length 

The first type of model protein considered here is the 
infinite straight cylinder originally introduced by West 
et al. [H]. (see Fig. [I] left), i.e., the protein cannot tilt. 
The axis of the cylinder is always aligned parallel to the 
z axis of the simulation box. Since this model protein 
spans the whole simulation box in the z direction and 
due to the periodic boundary conditions, the term "in- 
finitely long cylinder" is justified. The interaction of this 
simple model proteins and the lipid or solvent beads has 
repulsive contributions, which are described by a radially 
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shifted and truncated Lennard- Jones potential 

^''''^^ ' ~ 1 otherwise ' ^ > 



L [ot] 


Ept = 3 e 


Ept = 6 e 


4 


14.1 ± 2.1 


12.4 ± 2 


6 


12.5 ± 0.8 


11.9 ± 2.2 


8 


11.1 ± 1.2 


7.74 ± 1.15 



where r — + if^ denotes the distance of the interac- 
tion partners in the xy plane, a is given hy a = (at~\-Oi)l2 
for interactions with beads of type i (i = h, t, and s for 
head, tail, and solvent beads, respectively), CTq = ctj, and 
Vhj has been defined above (Eq. ([s])). The direct protein- 
protein interactions have the same potential (Eq. (4])) 
with a = at and ao — 2(Tf. In SI- units, the cylinders thus 
have a diameter of 2 nm, which roughly corresponds to 
the diameter of the 8-stranded /3-barrel OmpA [3T] . 

In addition, protein cylinders attract tail beads on a 
hydrophobic section of length L. This is described by 
an additional attractive potential that depends on the z 
distance between the tail bead and the protein center. 
The total potential reads 

Vptir, z) = ept {V,,p{r) + Kttr(r-) x M^p(z)) , (5) 

with the attractive Lennard- Jones contribution 

r ^ij(l)-14j(2) iir-ao<a 
Vattr(r) = <^ (^) - Vlj{2) if a < r - fxo < 2a , 
I otherwise 

(6) 

and a weight function Wp (z) , which is unity on a stretch 
of length 21 ~ L — 2at and crosses smoothly over to zero 
over a distance of approximately at at both sides. Specif- 
ically, we use 

1 if \z\ < I 

Wp = <! cos2 (§ (1^1 - /)) if ; < |z| < / + f . (7) 
otherwise 



TABLE I: Average tilt angle (q) of single spherocylinders in 
membranes. 



model bilayers. The distribution of solid angles assumes 
a maximum at tilt zero, and the average tilt angles (a) 
for the parameters used in the present work are less than 
15 degrees, as shown in Table [ij These tilt angles should 
be regarded as valid in the general case, since we are not 
aiming at modeling any specific protein. Our results are 
compatible with an earlier simulation study by Venturoli 
et aL[35], where the average tilt angle of proteins with 
both negative and positive mismatch of up to 40% was 
also only slightly affected by the mismatch. Only when 
a considerable positive mismatch of more than 70% was 
present, a significant increase in average tilt angle for 
model proteins of small aspect ratio (diameter to length) 
was observed. 

Since a protein is "rough" on an atomic scale, rep- 
resenting a complex structure like a protein as a sim- 
ple smooth cylindrical object may seem to be a rather 
crude approach. Nevertheless, our way of modeling the 
proteins can be justified by the fact that e.g. a-helices 
are packed with vanishing little free space within the he- 
lices [33] and are therefore fairly smooth on the scale of 
~ 10 A. There are no large cavities into which chains or 
even whole molecules would fit [34j. 



Simulation Method 



The hydrophobicity of the protein is tuned by the pa- 
rameter ept. 



2. Cylinder of finite length 

The second type of protein is also a straight cylinder, 
where the hydrophobic section of length L is modeled 
analogously to the infinite cylinder. But now the cylinder 
is allowed to tilt away from the z direction of the simula- 
tion box (Fig. [T] right) . The protein has finite length L 
and is capped at both ends by effectively repulsive ("hy- 
drophilic" ) hemispheres. More precisely, the protein is 
parametrized by a line of length L, and the interactions 
between lipids and proteins are given by Eq. ([s]) with 
"■q = ft as before, where r is now the shortest distance 
between the center of a lipid bead and the protein line. 
The interactions between two proteins are purely repul- 
sive and given by Q with ao = 2at and r the minimum 
distance between the two interacting protein lines. 

Even though these cylinders are allowed to tilt, they 
basically stay aligned with the membrane normal in our 



We have carried out Monte Carlo simulations at con- 
stant pressure P — 2.0 e/af and constant temperature T 
with periodic boundary conditions in a simulation box of 
variable shape and size. Thus, we are performing Monte- 
Carlo simulations in an NPT ensemble with effective 
Hamiltonian 



H^s^U + PV~ NkpT \n{V/Vo) 



(8) 



where U is the interaction energy, V the volume of the 
simulation box, Vq an arbitrary reference volume and N 
the total number of beads. 

In practice, three main types of Monte-Carlo moves 
were proposed and then accepted or rejected according 
to a Metropolis criterion, namely 1) translational local 
moves of the lipid beads, 2) protein translation and rota- 
tion moves, and 3) global moves which change the volume 
of the simulation box or its shape. Most of these moves 
have been discussed earlier [21] . Here we only sketch the 
algorithm for the new protein rotation moves. The new 
trial direction of the protein is generated in several 
steps. First a vector u, which lies randomly distributed 
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on a unit sphere, is generated. Choosing a random vec- 
tor on the surface of a unit sphere is efficiently done by 
applying an algorithm proposed by Marsaglia [33]: Two 
random variables ri and within the interval (—1,1) 



are generated and C — 



is calculated. If C > 1 



the random numbers are discarded and a new pair r\ , r2 
has to be generated. For < 1 the components of the 
random unit vector in Cartesian coordinates are given by 



(9) 



1 - 2C 



If the new direction of the protein was chosen to be u, 
the acceptance rate of the move would be prohibitively 
small. Therefore, u is scaled down by Atiit and the vector 
t = AtiitU is added to the tip of the normalized direction 
dp of the protein. Finally, normalizing the resulting vec- 
tor dp = dp + 1 to d'p gives the trial direction of the pro- 
tein used to calculate the energy change for the Metropo- 
lis criterion. By varying the scaling factor Atut during 
the "prerun" the acceptance rate of the tilt moves can 
be adjusted to the desired rate. Effectively, this means 
adjusting the opening angle of the cone around dp, which 
limits the maximum tilt angle Op^max of the move. 

During one Monte-Carlo step (MCS) there is on aver- 
age one attempt to move each bead and protein, and one 
attempt to rotate the protein. Since the global moves re- 
quire rescaling of all particle coordinates, which is rather 
expensive from a computational point of view, they are 
performed only every 50th MCS on average. 



D. Elastic Theory 

In the following, we briefly sketch the elastic theory 
with which wc compare the simulation results [l0 l [T4) , [T8]. 
and which has proved to describe very well the proper- 
ties of our model membranes in the fluid state [IH] ■ The 
membrane is treated as a system of coupled monolayers, 
which fluctuate subject to the constraint that the vol- 
ume of lipids is locally conserved. It is also considered 
to be almost flat, i.e., the positions of both monolayers 
can be parametrized by single- valued functions hi(x,y). 
Here we are mainly interested in the thickness deforma- 
tions, which we characterize by the locally smoothed de- 
formation profiles $ci(a^, y) of monolayers about the un- 
perturbed value Iq ( i.e., the total monolayer thickness is 
io + ^'ci)- The free energy of monolayer thickness defor- 
mations is then written as il5i 



F 
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all 
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all 


0.15 at'^ 


kA/tl 


all 


1.3 e/ai 
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all 


-0.26 e 


tR 


4 (Tt 


-0.94 at 


tR 


6 (Tt 


0.3 at 


tR 


8 at 


1.44 at 


Co 


4 at 


-0.11 


Co 


6 at 


0.05 a^^ 


Co 


8 at 


0.22 at^ 



TABLE II: Elastic parameters used in the theoretical calcu- 
lations. Taken from Ref. 1181. 



where and are the bending and compressibility 
modulus of the bilayer, cq is the spontaneous curvature, 
C a parameter related to the spontaneous curvature, and 
kc the Gaussian curvature. All these parameters have 
been determined for our model membranes in earlier work 
from analyses of the membrane fluctuations and stress 
profiles ^8j. We note that the terms cq and kc can 
be rewritten as pure surface terms. Protein inclusions 
distort the membranes by imposing a fixed thickness at 
their boundaries [TU], <I>ei = tn, and may introduce ad- 
ditional surface fields, which enter in the same way as 
the spontaneous curvature term [181 and thus effectively 
renormalize cq . Hence the effect of proteins on the mem- 
brane can be characterized by the two parameters Ir and 
Co, the renormalized curvature. They have been deter- 
mined for our model in earlier work |18j from the dis- 
tortion profiles of single membranes about infinitely long 
straight cylinders. Table [TT| summarizes the values of all 
elastic parameters for the simulation models studied in 
this work. 



^(A<i>e 



A:Gdet(%$oi) \ , 



The elastic free energy, Eq. ( 10 ), can be minimized an- 
alytically for membranes containing only a single cylin- 
drical protein. For two proteins, it must be minimized 
numerically. To this end, the free energy integral was 
discretized in real space using a square grid of spatial 
discretization h = 0.25 (Xt, a system size of 50 x 50 and 
a second-order difference scheme to evaluate the deriva- 
tives. The boundary condition was implemented by set- 
ting (j) = tji inside the inclusion. The minimization was 
done via a steepest descent method, using a relaxation 
scheme [36'. This procedure gave deformation profiles 
with which we could compare the simulation results (see 
below) . 



III. SIMULATION RESULTS 

A. Lipid-mediated interactions between inclusions 

We now turn to discussing our simulation results for 
the membrane-induced interactions between cylindrical 
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inclusions in the bilayer. In this subsection, we focus 
on the effect of protein orientation fluctuations on the 
effective interactions between the model proteins, i.e., 
the PMF. 

The radial distribution function g(r) as a function of 
the protein-protein distance r was obtained from sim- 
ulation runs using the technique of successive umbrella 
sampling [37 combined with a reweighting procedure. As 
starting configurations we used equilibrated systems with 
750 to 760 lipids and two simple transmembrane proteins 
of diameter Sat. A first estimate of g{r) was obtained 
during pre- runs of length 2 x 10^ MCS. Then, biased runs 
of 3 X 10^ MCS were performed to improve the statis- 
tics of configurationally less frequent protein-protein dis- 
tances. In the case of tiltable proteins, 10 x 10® MCS had 
to be performed to achieve comparable accuracy in the 
histogram data. 

After removing the bias from these results and com- 
bining the overlapping distributions the effective poten- 
tial w^r) — ^kBT\ng{r) can be extracted. At small 
inclusion-inclusion distances direct interactions of pro- 
teins and, in the case of long proteins, additional de- 
pletion induced attraction due to the solvent particles 
come into play. Since our main interest lies in the study 
of lipid-mediated medium and long ranged interactions, 
these parts of the curves have been cut off. Our sam- 
pling procedure resulted in a statistical error of about 
±0.3e within each umbrella window. The distance range 
between 4 and 10 at was typically covered by 7-8 windows 
(with larger windows at larger distances) , which gives an 
accumulated error of ±2 e in the first minimum of the 
PMFs in the worst case, and of at most ±1 e at distance 
6 a. 

The hydrophobicity of the protein, i.e., the interaction 
strength between the hydrophobic core of the membrane 
and the hydrophobic part of the inclusion, is tuned by 
the parameter Cpt . Only for sufficiently high values of Cpt 
do the model proteins locally distort the bilayers. In the 
following we will consider two cases, weakly hydropho- 
bic proteins with ept = 3.0 e and strongly hydrophobic 
proteins with ept = 6.0 e. In the first case lipid-protein 
contact is approximately energetically equivalent to lipid- 
lipid contact, whereas in the second case lipid-protein 
contact is highly preferred. Weakly hydrophobic pro- 
teins only induce small perturbations of the lipid envi- 
ronment, strongly hydrophobic proteins lead to strong 
deformations. 

Furthermore, the hydrophobic length L has been var- 
ied from L — Aat to L = 8 at. The hydrophobic thickness 
of the membrane is 6 at, hence L — Aat corresponds to 
a negatively mismatched protein, and L ~ 8at to a pos- 
itively mismatched protein. 

The simulation results are compiled in Figs. [2j [3j and 
|4j The left and right graphs in each figure correspond to 
straight and tiltable cylinders, respectively. One immedi- 
ately notices the first effect of orientational fluctuations: 
Due to lipid packing in the vicinity of the inclusions the 
curves feature oscillations with a wavelength of approx- 
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FIG. 2: Potential of mean force of two isolated proteins 
as a function of the protein distance at medium (solid) and 
strong (dashed) hydrophobicity ept for proteins with negative 
hydrophobic mismatch {L = 4 at). The proteins of the left 
graph were always aligned parallel to the z axis of the simu- 
lation box, whereas proteins in the right graph were allowed 
to tilt. 
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FIG. 3: Potential of mean force of two isolated proteins as a 
function of the protein distance at medium (solid) and strong 
(dashed) hydrophobicity ept for hydrophobically matching 
proteins (L — 6at). (left: aligned protein, right: tiltable 
proteins). 
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FIG. 4: Potential of mean force of two isolated proteins 
as a function of the protein distance at medium (solid) and 
strong (dashed) hydrophobicity ept for proteins with positive 
hydrophobic mismatch (L — 8 at) (left: aligned protein, right: 
tiltable proteins). 
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imately I at, i.e. the diameter of the hpid tail beads. 
The higher ept the more pronounced the hpid layering. 
This holds both for straight and tiltable proteins. How- 
ever, the oscillations are significantly reduced in the lat- 
ter case. Hence orientation fluctuations reduce the effect 
of lipid layering on the protein-protein interactions. 

Next we discuss separately the PMFs for strongly and 
weakly hydrophobic proteins. 



1. Strongly hydrophobic proteins 

A general feature that straight and tiltable protein 
have in common at Cpt = 6.0 e is an enhanced tendency 
to aggregate for hydrophobically mismatched proteins, 
both for positive or negative mismatch (dashed curves in 
Figs. [2] and |4|. The effective interactions between mis- 
matched proteins have an attractive long-range contribu- 
tion which is not detectable for hydrophobically matching 
proteins (Fig. [S]). Compared to straight proteins, how- 
ever, the range of this additional attraction is reduced 
for proteins with orientational fluctuations, and the at- 
traction essentially vanishes within the error at distances 
beyond r ~ 7 ctj . 

The first and deepest attractive minimum of the PMFs 
slightly below r = 4at lies at approximately the same 
energy for both types of models, i.e. between —17 and 
— 15e for proteins of length L = Aat and 6 at and about 
—40 e for proteins of length L — Sat- This minimum 
is related to lipid bridging and a consequence of the 
strong attraction between lipids and proteins. At dis- 
tance r ^ Aat, one layer of lipids is in contact with both 
proteins, thus stabilizing the conformation at that dis- 
tance. The lipid bridging energy increases with increas- 
ing hydrophobic contact area between lipid tails and pro- 
teins. Compared to hydrophobically matching proteins, 
it is thus reduced for negatively mismatched proteins 
and enhanced for positively mismatched proteins. In ad- 
dition, both negatively and positively mismatched pro- 
teins have the above-mentioned attractive interactions 
due to the elastic deformation of the layer. Compared 
to hydrophobically matching proteins, the contributions 
of lipid bridging and hydrophobic mismatch are oppo- 
site for negatively mismatched proteins, and cumulative 
for positively mismatched proteins. As a result the first 
minimum has approximately the same energy for pro- 
teins with negative hydrophobic mismatch and matching 
proteins, whereas proteins with positive mismatch show 
the strongest attraction. 

2. Medium hydrophobicity 

For weakly hydrophobic proteins with ept — 3.0 e, the 
effect of hydrophobic mismatch at short distances is much 
less pronounced than in the case of strongly hydrophobic 
proteins. This is especially true for proteins with fluctu- 
ating orientations. In the latter case, the first minimum 



has almost the same depth for all types of mismatch, and 
the PMFs for different L are not further apart than 3 e 
(solid curves in the right graphs of Fig. [5] [Sj and|4|. 

At intermediate distances, however, the interaction be- 
tween proteins with fluctuating orientations is found to 
depend crucially on the type of mismatch: For positively 
mismatched proteins, it has a long range attractive con- 
tribution similar to that observed for strongly hydropho- 
bic proteins. For negatively mismatched proteins, this 
contribution is still present at distances r < 7 at, but 
then it turns around and gives way to a repulsive in- 
teraction. The repulsive barrier is small, but significant 
within the error of ±0.5 e at distance 8 at- The behav- 
ior for matching proteins is intermediate - the layering 
interactions are superimposed by an attractive contribu- 
tion up to r ~ 6tTt, which disappears within the error at 
larger distances. We note that this behavior is in marked 
contrast both to the behavior of strongly hydrophobic 
proteins and of weakly hydrophobic, but straight inclu- 
sions. In the latter cases, one observes a purely attrac- 
tive long-range hydrophobic mismatch interaction, which 
vanishes for hydrophobically matched proteins. 

In fact, the elastic theory does predict a repulsive bar- 
rier in the interaction free energy of two transmembrane 
proteins for positive and negative hydrophobic mismatch 
[TFj . The maximum of this small repulsive barrier should 
be found in the range of r 6 at to 8 at and can be 
associated with a peak, i.e., a soft mode, in the spec- 
trum of thickness fluctuations of the pure bilayer. One 
may speculate that the repulsive shoulders observed here 
are caused by the same effect. However, our simulations 
results show no repulsive barrier for proteins of length 
L — 8.0 at. The attractive range for positively mis- 
matched proteins rather extends up to a protein-protein 
separation of about 9 at and levels off without any observ- 
able repulsive region. Furthermore, the attractive behav- 
ior of hydrophobically matching proteins is not predicted 
by the elastic theory either. 

Other comparable simulations on tiltable proteins have 
been performed by Schmidt et al. Both for negative 
and positive hydrophobic mismatch they observe a long- 
range, lipid-mediated attraction. In the case of zero mis- 
match this interaction is less attractive [35] ■ Further, 
their PMFs exhibit an oscillating " flne structure" , which 
they attribute to the discreteness of the membrane. 

In another study, de Meyer et al. have compared 
PMFs for proteins which were allowed to tilt, and pro- 
teins which were always aligned parallel to the z-axis [17] . 
Their results for proteins comparable to our capped cylin- 
drical inclusions of diameter Dp = 3at can be summa- 
rized as follows: In the case of positive mismatch the 
long-range interaction between two proteins should be 
influenced by the degree of protein tilt. If the protein 
was not allowed to tilt, they found an attractive region 
for small distances, followed by a small repulsive barrier 
and a shallow and broad minimum at large protein sep- 
aration. If tilting of proteins was allowed, only a small, 
but broad repulsive barrier at medium inter-protein dis- 
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tances was found after the typical attractive range at 
small distances for both positive and negative mismatch. 
The shallow attractive region at larger distance has van- 
ished. In our model of long cylindrical proteins without 
tilt a shallow repulsive barrier at medium distances was 
also present, but only at very weak hydrophobicity of 
ept = 1.0 e [T5]. Other effects due to lipid layering were 
much less pronounced in the study of de Meyer et al., 
since they were using soft, merely repulsive potentials. 
But at very small protein separations they did observe 
slight oscillations in the PMFs, which they assigned to 
the free energy needed to remove lipids in between the 
proteins. 



B. Profiles around two proteins 

After having discussed the effective lipid-mediated in- 
teractions between two transmembrane inclusions, we 
now turn to the closer investigation of structural lipid 
rearrangements around two model proteins at small and 
medium separation. In this section, we will consider 
strongly hydrophobic proteins with ept = 6 e. 

As we have seen in the previous section, one character- 
istics of our lipid model with its hard-core, Lennard- Jones 
type interactions is the presence of lipid packing phenom- 
ena, which are usually much less prominent in systems 
modeled with soft dissipative particle dynamics poten- 
tials |18j . Especially in the area surrounding proteins 
one can expect ordering of the lipids, which is rather dis- 
tinct from the situation further away from inclusions in 
the lipid bulk, and which may affect the interaction with 
other inclusions. Here, two-dimensional profiles give use- 
ful insight into the structure of membrane-characterizing 
quantities. To determine them, we have simulated equi- 
librated bilayers containing two capped proteins, keeping 
the positions of the proteins fixed, but allowing for tilt 
fluctuations. As before, we compare systems with cylin- 
drical proteins of varying hydrophobicity, hydrophobic 
mismatch and protein-protein distance. 

The two-dimensional profiles can be separated into two 
regions: In the inner region between the two inclusions, 
the profiles reflect the combined effect of both proteins. 
In the outer region, the influence of one protein is ge- 
ometrically screened, and the profiles mainly reproduce 
the perturbations induced by one single protein. 

1. Areal tail bead density 

First, the areal density of tail beads Xt around two pro- 
teins is investigated. In the following graphs, the values 
are always given with respect to one monolayer, i.e. the 
average number of tail beads obtained for each bin lying 
in the xy plane was divided by 2. The reference value 
of Xt in the unperturbed bilayer is simply the number of 
tail beads per lipid divided by the average area per lipid 
at r = 1.3 ksT, i.e. Xt = 6/1.38(7^2 = 4.3 crr^- 



d = 3.99 at d = 6.45 at 




20 24 2B 32 20 24 2B 32 



FIG. 5: Areal tail density profiles for systems containing two 
translationally immobile hydrophobically matched proteins 
(L — 6 at) for two different protein distances, obtained by 
binning within a quadratic grid with bin size 0.5 x 0.5 af. 
The profiles reflect the underlying grid. 



Both length and height of a single bin is set to 0.5 at. 
Therefore each sampled cell has an area of 0.25 only. 
At this resolution interesting details of the lipid structure 
in the vicinity of the inclusions are revealed. 

Not surprisingly, negative hydrophobic mismatch, 
which goes along with thinning of the surrounding bi- 
layer, leads to a reduction of the areal tail density, 
whereas positive hydrophobic mismatch enhances the 
areal tail density along with the membrane thickness. 
The fine structure exhibits additional features. In the 
following, we will focus on hydrophobically matching pro- 
teins of length L = 6 CTt . The layering of the lipids around 
the proteins, which gave rise to the oscillatory behavior of 
the PMFs in the previous section, manifests itself clearly 
in concentric rings with enhanced (light) and reduced 
(dark) areal tail bead density (Fig. [5]). The variations in 
areal density are usually in the order of 20% compared to 
the average value at larger distance from the inclusions. 
If the proteins are close to each other (Fig. [5] left), they 
pin the lipids entirely i.e., the areal tail density exhibits 
sharp local peaks at two well-defined positions between 
the two proteins. 



2. Thickness profiles 

Next, we inspect the thickness profiles around two 
cylindrical inclusions at fixed distances. Specifically, we 
consider short distances, slightly below a protein-protein 
separation of ^ 4 dt , medium distances at a separation of 
~ 6.5 at, and larger separations lying between 8 and 9 at. 
A regular grid of cells with an area of Ax Ay — 1.0 cr^ 
was used to subdivide the system. Averaging of the 
membrane thickness was done within each quadratic bin 
thus obtained. If necessary, the systems were rotated 
such that the connecting line between the centers of the 
proteins was parallel to the x-axis for better comparison 
and visualization purposes. Light shading indicates thick 
membranes, dark shading thin membranes. For compar- 
ison, we have also calculated the same profiles with the 
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FIG. 6: Thickness profiles around two translationally immo- 
bile cylindrical proteins with negative hydrophobic mismatch 
(L — 4.0 at). Left: Simulations, Right: Elastic theory 



FIG. 7: Thickness profiles around two translationally immo- 
bile cylindrical proteins with positive hydrophobic mismatch 
{L — 8.0 at). Left: Simulations, Right: Elastic theory 



elastic model described earlier (Sec. |IID ), using the elas- 
tic parameters obtained in Ref. 18 (Table The re- 
sults are shown in Figs. [6j [Tj andjSj 

Both for small and medium protein-protein separations 
the perturbations induced by the proteins overlap. This 
amplifies the reduction or increase of thickness which 
would be caused by a single protein of negative or pos- 
itive hydrophobic mismatch, respectively. At large dis- 
tances, the proteins are surrounded by separate pertur- 
bation shells. 

Fig. [6] shows the results for proteins with negative mis- 
match. For small and medium protein-protein distances 
the bilayer thickness is smallest in the region between 
the proteins (Fig. [6| top left and middle left), which is in 
remarkable agreement with the predictions of the elastic 
theory (Figs. [6] top right and middle right). For protein- 
protein distances 5 — 7<7t i the elastic theory even predicts 
a certain amount of " undershooting" , i.e., the membrane 
thickness on the axis between the two proteins is pre- 
dicted to be smaller than at the protein surface. This is 
also found in simulations, albeit to a much lesser extent. 

Proteins with positive hydrophobic mismatch {L = 
8.0 at) are surrounded by a clearly thickened bilayer 
(Fig. [t]). As for negatively mismatched proteins, the 
membrane perturbations are enhanced in the region be- 
tween the proteins. For protein distances 5 — 7 at, the 



elastic theory predicts "overshooting", i.e., the thick- 
ness is expected to be largest at the point right between 
the two proteins. However, the simulation data do not 
show this effect. In this respect, positively mismatched 
proteins show a different behavior than negatively mis- 
matched proteins. 

Finally, the bilayer thickness between hydrophobically 
matching proteins {L = 6.0 at) remains practically unal- 
tered as expected, except in the very close vicinity of the 
proteins, where a small annulus of enhanced thickness 
develops (Fig. [s] left). Elastic theory predicts a slightly 
enhanced thickness in the vicinity of the proteins (Fig. |8j 
right). 



IV. SUMMARY AND DISCUSSION 

In summary, we have studied protein-protein interac- 
tions in model lipid membranes, using two types of model 
proteins: Infinitely long cylinders with fixed orientation, 
and finite spherocylinders with fiuctuating orientation. 
In general, the interactions are dominated by packing ef- 
fects and smooth attractive hydrophobic mismatch inter- 
actions. Besides the geometric mismatch the hydropho- 
bicity of the proteins is a key quantity. 

The comparison of the PMFs for the two protein types 
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FIG. 8: Thickness profiles around two translationally immo- 
bile cylindrical proteins without mismatch (L — 6.0 at). Left: 
Simulations, Right: Elastic theory. 

allowed to assess the effect of orientation fluctuations, ft 
can be summarized as follows: The main and most promi- 
nent effect of orientation fluctuations is to reduce the 
lipid packing effects on the PMFs, i.e., the correspond- 
ing short range oscillatory behavior. Furthermore, they 
slightly reduce the range of the hydrophobic mismatch 
interactions. For weakly hydrophobic proteins, they may 
also affect qualitatively the shape of the potentials: fn 
the case of negatively mismatched proteins, orientation 
fluctuations introduce a weak repulsive maximum in the 
PMFs which is not present for orientationally fixed pro- 
teins. This is consistent with earlier results by de Meyer 
et al. jl7| , using a different protein model with hydropho- 
bic interactions that would be considered 'weak' in our 
context. 

The direct experimental measure of interaction poten- 
tials between membrane proteins is still a challenging 
task due to the required temporal and spatial resolution. 
The "protein" diameter in our study roughly corresponds 
to that of gramicidin, for which the effect of hydropho- 
bic mismatch interactions was verified experimentally in 
1999 by Harroun and coworkers [35] • The hydrophobic 
mismatch effect has also been investigated by system- 
atic studies of synthetic model peptides [2j [40] . In these 
studies, the evidence for mismatch-induced interactions 
was mainly based on the observation of protein cluster- 
ing. Recently, spatially resolved interaction energies be- 



tween proteins of positive hydrophobic mismatch have 
been calculated for mobile ATP-synthase c-rings in mem- 
branes basically consisting of phosphatidyl glycerophos- 
phate from the probability distribution of the center- 
to-center distance by Casuso et al. j4lj. After a soft- 
core short range repulsion attributed to the structural 
perturbation of the lipid/protein organization, a mini- 
mum could be observed in the interaction energy. The 
authors consider this attractive potential to be caused 
by membrane-mediated interactions due to the deforma- 
tions of the lipid bilayer. The ATP-synthase c-rings of 
diameter 65 A ± 5 A would correspond to proteins of di- 
ameter 10 at in our systems. Therefore, direct compari- 
son with our simulation data is not possible. Neverthe- 
less, the experimental finding of protein-protein attrac- 
tion by mismatch-induced deformation is in accordance 
with the resulting smooth attractive interaction in our 
simulations. 

For the case of strong interactions, we have also studied 
the distribution of lipid tails in systems containing two 
proteins and found that two proteins can form complexes 
with surrounding lipids, pinning these to their surface. 
We speculate that such effects may also occur in biolog- 
ical environments. Furthermore, we have compared the 
thickness profiles in systems containing two proteins with 
the elastic theory of two coupled monolayers. We found 
remarkable agreement between the theoretical prediction 
and the simulation results, with only slight qualitative 
deviations at intermediate protein distances ^ 7 at- here, 
the elastic theory predicts "over"- or "undershooting" ef- 
fects which are much weaker or not discernible at all in 
the simulations. The theoretical prediction for the PMFs 
has not been shown here (see [HI), because it is difficult 
to compare with the simulation data due to the lipid 
packing effects. Regarding the smooth hydrophobic mis- 
match part, the theory is in rough qualitative agreement 
with the simulations, except that its predicts a repulsive 
barrier at distances ~ 8 (Tt which is not observed in the 
simulations. The discrepancy is possibly related to the 
slight deviations between theory and simulations in the 
thickness profiles mentioned above. 



Acknowledgement 

Provision of computing resources by the HLRS 
(Stuttgart), NIC (Jiihch), and PC2 (Paderborn) is grate- 
fully acknowledged. The configurational snapshots were 
visualized using VMD . This work was funded in part 
by the DFG within the Sonderforschungsbereich SFB 
613, SFB 625, and SFB TR 6. 



[1] Berg, J. M.; Tymoczko, J. L.; Stryer, L. Biochemistry, 
5th ed.; W. H. Freeman and Company, 2002. 

[2] de Planque, M. R. R.; Killian, J. A. Protein-lipid inter- 
actions studied with designed transmembrane peptides: 



role of hydrophobic matching and interfacial anchoring 
(Review). Molecular Membrane Biology 2003, 20, 271- 
284. 

[3] May, S. Theories on structural perturbations of lipid bi- 



10 



layers. Current Opinion in Colloid & Interface Science [24 
2000, 5, 244-249. 

Coarse- Graining of Condensed Phase and Biornolecular 
Systems; Voth, G. A., Ed.; CRC Press, Boca Raton, [25 
2009. 

Miiller, M.; Katsov, K.; Schick, M. Biological and syn- 
thetic membranes: What can be learned from a coarse- [26 
grained description? Phys. Rep. 2006, 434, 113-176. 
Descrno, M. Mesoscopic Membrane Physics: Concepts, 
Simulations, and Selected Applications. Macromolecular 
Rapid Communications 2009, 30, 752-771. 
Schmid, F. Toy amphiphiles on the computer: What can [27 
we learn from generic models? Macromolecular Rapid 
Commumcations 2009, 30, 741-751. 

Marcolja, S. Lipid-mediated protein interactions in mem- [28 
branes. Biochimica et Biophysica Acta 1976, 4-55, 1-7. 
Dan, N.; Pincus, P.; Safran, S. A. Membrane-Induced In- 
teractions between Inclusions. Langmuir 1993, 9, 2768- [29 
2771. 

Aranda-Espinoza, H.; Berman, A.; Dan, N.; Pincus, P.; 
Safran, S. Interaction between inclusions embedded in [30 
membranes. Biophysical Journal 1996, 71, 648. 
May, S.; Ben-Shaul, A. A molecular model for lipid- 
mediated interaction between proteins in membranes. [31 
Physical Chemistry Chemical Physics 2000, 2, 4494- 
4502. 

Lagiic, P.; Zuckermarm, M. J.; B., R. Lipid-Mcdiated In- [32 
teractions between Intrinsic Membrane Proteins: A The- 
oretical Study Based on Integral Equations. Biophysical 
Journal 2000, 79, 2867. 

Bohninc, K.; Kralj-Iglic, V.; May, S. Interaction between 

two cylindrical inclusions in a symmetric lipid bilayer. [33 
Journal of Chemical Physics 2003, 119, 7435-7444. 
Brannigan, G.; Brown, F. L. H. A Consistent Model [34 
for Thermal Fluctuations and Protein-Induced Deforma- 
tions in Lipid Bilayers. Biophysical Journal 2006, 90, 
1501. [35 
Brannigan, G.; Brown, F. L. H. Contributions of Gaus- 
sian Curvature and Nonconstant Lipid Volume to Pro- 
tein Deformation of Lipid Bilayers. Biophysical Journal [36 
2007, 92, 864 876. 

Schmidt, U.; Guigas, G.; Weiss, M. Cluster Formation [37 
of Transmembrane Proteins Due to Hydrophobic Mis- 
matching. Physical Review Letters 2008, 101, 128104. 
de Meyer, F. ,J.-M.: Vcnturoli, M.; Smit, B. Molecular [38 
Simulations of Lipid-Mediated Protein-Protein Interac- 
tions. Biophysical Journal 2008, 95, 1851-1865. [39 
West, B.; Brown, F. L. H.; Schmid, F. Membrane-Protein 
Interactions in a Generic Coarse-Grained Model for Lipid 
Bilayers. Biophysical Journal 2009, 96, 101-115. 
de Meyer, F. J.-M.; Rodgcrs, J. M.; Willems, T. F.; 
Smit, B. Molecular Simulation of the Effect of Choles- [40 
terol on Lipid-Mcdiated Protein-Protein Interactions. 
Biophysical Journal 2010, 99, 3629-3638. 
Goulian, M. Inclusions in membranes. Current Opinion 
in Colloid and Interface Science 1996, 1, 358-361. 
Weikl, T. R. Fluctuation-induced aggregation of rigid [41 
membrane inclusions. Europhysics Letters 2001, 54, 547. 
Reynwar, B. J.; lUya, G.; Harmandaris, V. A.; 
Miiller, M. M.; Kremer, K.; Deserno, M. Aggregation and [42 
vesiculation of membrane proteins by curvature-mediated 
interactions. Nature 2007, 447, 461. 
[23] Lenz, O.; Schmid, F. A simple computer model for liquid 
lipid bilayers. J. Mol. Liquids 2005, 117, 147-152. 



Schmid, F.; Diichs, D.; Lonz., O.; West, B. A generic 
model for lipid monolayers, bilayers, and membranes. 
Computer Physics CornmumcaMons 2007, 177, 168-171. 
Lenz, O.; Schmid, F. Structure of Symmetric and Asym- 
metric "Ripple" Phases in Lipid Bilayers. Physical Re- 
view Letters 2007, 98, 058104. 

West, B.; Schmid, F. In NIC Symposium 2010; 
Miinstcr, G., Wolf, D., Kremer, M., Eds.; Forschungszon- 
trum Jiilich, 2010; Chapter Membrane-Protein Inter- 
actions in Lipid Bilayers: Molecular Simulation versus 
Elastic Theory, pp 279-286. 

Neder, J.; West, B.; Nielaba, P.; Schmid, F. Coarse- 
Grained Simulations of Membranes under Tension. Jour- 
nal of Chemical Physics 2010, 132, 115101. 
Langel, U.; Cravatt, B. F.; Graslund, A.; von Heijne, G.; 
Land, T.; Niesscn, S.; Zorko, M. Introduction to Peptides 
and Proteins; CRC Press, 2010. 

Wimley, W. C; White, S. H. Experimentally determined 
hydrophobicity scale for proteins at membrane interfaces. 
Nature Structural Biology 1996, 3, 842-848. 
Bechingcr, B. Towards Membrane Protein Design: pH- 
sensitive Topology of Histidiuc-containing Polypeptides. 
Journal of Molecular Biology 1996, 263, 768-775. 
Alberts, B.; Johnson, A.; Lewis, J.; Raff, M.; Roberts, K.; 
Walter, P. Molecular Biology of the Cell; Garland Sci- 
ence, 2002. 

Vcnturoli, M.; Smit, B.; Sporotto, M. M. Simulation 
Studies of Protein-Induced Bilayer Deformations, and 
Lipid-Induced Protein Tilting, on a Mesoscopic Model 
for Lipid Bilayers with Embedded Proteins. Biophysical 
Journal 2005, 88, 1778. 

Handbook of Molecular Biophysics; Bohr, H. G., Ed.; 
Wiley- VCH Vorlag GmbH & Co. KGaA, 2009. 
Scott, H. L.; Coc, T. ,J. A Theoretical Study of 
Lipid-Protein Interations in Bilayers. Biophysical Jour- 
nal 1983, 42, 219-224. 

Marsaglia, G. Choosing a point from the surface of a 
sphere. Annals of mathematical statistics 1972, 43, 645- 
646. 

Schmid, F.; Schick, M. Liquid phases of Langmuir mono- 
layers. Journal of Chemistry Physics 1995, 102, 2080. 
Virnau, P.; Miiller, M. Calculation of free energy 
through successive umbrella sampling. Journal of Chem- 
ical Physics 2004, 120, 10925-10930. 
Schmidt, U.; Guigas, G.; Weiss, M. Schmidt, Guigas and 
Weiss Reply. Physical Review Letters 2009, 102, 219802. 
Harroun, T. A.; Heller, W. T.; Weiss, T. M.; Yang, L.; 
Huang, H. W. Experimental Evidence for Hydropho- 
bic Matching and Membrane-Mediated Interactions in 
Lipid Bilayers Containing Gramicidin. Biophysical Jour- 
nal 1999, 76, 937-945. 

Sharpe, S.; Barber, K. R.; Grant, C. W. M.; 
Goodyear, D.; Marrow, M. R. Organization of model he- 
lical peptides in lipid bilayers: Insight into the behaviour 
of single-span protein transmembrane domains. Biophys- 
ical Journal 2002, 83, 345-358. 

Casuso, I.; Sens, P.; Rico, F.; Scheuring, S. Experimental 
Evidence for Membrane-Mediated Protein-Protein Inter- 
action. Biophysical Journal 2010, 99, L47-L49. 
Humphrey W.; Dalke, A.; Schultcn, K. VMD Vi- 
sual Molecular Dynamics. Journal of Molecular Graphics 
1996, 14, 33-38. 



